Astronomy Sz Astrophysics manuscript no. absAGN'last 


© ESQ 2008 


February 4, 2008 





An X-ray view of absorbed INTEGRAL AGN 

A. De Rosa 1 , L. Bassani 2 , P. Ubertini 1 , F. Panessa 1 , A. Malizia 2 , A. J. Dean 3 , and R. Walter 4 

1 INAF/IASF, via del Fosso del Cavaliere 100, F00133 Rome, Italy 

2 INAF/IASF, via Gobetti 101,1-40129 Bologna, Italy 

3 School of Physics and Astronomy, University of Southampton Highfield, Southampton, United Kingdom, 

4 INTEGRAL Science Data Centre, Chemin d'Ecogia 16, CH-1290 VERSOIX, Switzerland 

Received /Accepted 

ABSTRACT 

Aims. We present a 0.2-200 keV broad-band study of absorbed AGN observed with INTEGRAL, XMM-Newton, Chandra 
and ASCA to investigate the continuum shape and the absorbing/reflecting medium properties. 

Methods. The sources are selected in the INTEGRAL AGN sample to have a 20-100 keV flux below 8xl0~ n 
erg cm -2 s _1 (5 mCrab), and are characterized by a 2-10 keV flux in the range (0.8-10) x 10 -11 erg cm -2 s _1 . The 
good statistics allow us a detailed study of the intrinsic and reflected continuum components. In particular, the analysis 
performed on the combined broad-band spectra allow us to investigate the presence of Compton reflection features and 
high energy cut-off in these objects. 

Results. The column density of the absorbing gas establishes the Compton thin nature for three sources in which a 
measure of the absorption was still missing. The Compton thin nature of all the sources in this small sample is also 
confirmed by the diagnostic ratios F x /F[OIII]. The Compton reflection components we measure, reflection continuum 
and iron line, are not immediately compatible with a scenario in which the absorbing and reflecting media are one and 
the same, i.e. the obscuring torus. A possible solution is that the absorption is more effective than reflection, e.g. under 
the hypothesis that the absorbing/reflecting medium is not uniform, like a clumpy torus, or that the source is observed 
through a torus with a very shallow opening angle. The high energy cut-off (a lower limit in two cases) is found in 
all sources of our sample and the range of values is in good agreement with that found in type 1 Seyfert galaxies. At 
lower energies there is clear evidence of a soft component (reproduced with a thermal and/or scattering model), in six 
objects. 
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ing/absorbing medium are one and the same (the putative 
"molecular torus"). A cold iron line at 6.4 keV is often as- 
sociated to this continuum Compton reflection hump. 

The IBIS gamma-ray imager on board INTEGRAL 
is surveying the whole sky above 20 keV with a mCrab 
(~10 -11 erg cm -2 s _1 ) sensitivity in well exposed re- 
gions. With its angular resolution of 12 arcmin and the 
point source location accuracy of 1-2 arcmin for moder- 
ately bright sources (Ubertini et al. 2003) INTEGRAL has 
already detected a number of AGN, a large fraction of 
which were previously unknown as high energy emitters 
(jBird et al. 20061 [Bassani et al. 2006allBird et al. 2007|) . 

In this paper we present a deep broad-band spectral 
analysis of seven absorbed AGN detected by INTEGRAL. 
The sample is extracted from the Bassani et. al (2006a) sur- 
vey, updated to include a number of optical classifications 
obtained afterwards (see Bassani et al. 2006bp . The sample 
includes Seyfert 2 galaxies having a 20-100 keV flux be- 
low 5 mCrab and with X-ray data available at the time 
that this work started; from the sample we have excluded 
sources already studied by Beppo-SAX except for one ob- 
ject (ESO 103-G35) which was retained in order to provide 
a direct comparison between this work and previous stud- 
ies using Beppo-SAX data and therefore an a posteriori 
check of our analysis procedure which is affected by the 
limitation of using non simultaneous X and soft gamma- 
rays data. Overall we can conclude that the sample used 



1. Introduction 

The broad-band properties of Seyfert 2 AGN have been 
extensively studied in the past years mainly using Beppo- 
SAX data (IRisaliti 2002p . The continuum shape in these 
objects can be reproduced by a steep power-law with 
photon index T ~ 2. This model matches well a 
Comptonization scenario in which seed photons from a 
cold gas (the accretion disc) a re up-scattered in a hot 
corona of rclativistic electrons (Haardt & Maraschi 1991, 



Haardt, Maraschi & Ghisellini 1997). A photoelectric cut- 
off is present at X-ray energies and its value depends on 
the absorbing column density of the source. For N# in the 
10 22 — 10 23 cm -2 range the cut-off is below 2 keV, for N# in 
the 10 23 — 10 24 cm -2 range the cut-off is in 3-10 keV band 
and for N# > a^} — 1.5 x 10 24 cm~ 2 (i.e. the Compton 
thick regime with er^ 1 being the inverse of the Thomson 
cross section), the continuum can be observed only through 
its reflected component above 10 keV. 

The presence of a high energy cut-off is still an open 
issue. Risaliti (2002) found that the high energy cut-off is 
present only in ~ 30 per cent of a sample of 20 bright Sy 
2 observed with Beppo-SAX. Above 10 keV the presence of 
a Compton reflection component has often been observed, 
however is its origin still unclear as well as if the reflect- 
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in this study is representative of the population of type 2 
AGN detected by INTEGRAL above 10 keV; further broad 
band X/gamma-ray analysis on a complete sample of type 
2 AGN selected in the 20-40 keV band is on going and 
will be the objective of a future papeiQ. Four sources (IGR 
J12391-1610=LEDA170194 hereafter, IGR J07565-4139, 
IGR J12026-5349 and IGR J10404-4625) have been discov- 
ered first in soft gamma-ray band and immediately after ob- 
served with Chandra (LEDA170194, IGR J07565-4139 and 
IGR J12026-5349) and XMM-Newton (IGR J10404-4625). 
Two objects (ESO 103-G35 and IC 4518A) are known AGN 
(although IC 4518A was never reported at X-ray ener- 
gies before) and follow up observations with XMM-Newton 
have been used in our broad-band analysis. For one ob- 
ject (NGC 788) archival ASCA data have been used in 
this work. Using X-ray measurements in conjunction with 
INTEGRAL data allow us to build a spectrum in three 
decades of energy, and to distinguish between the various 
spectral components characterizing the broad-band spectra 
of absorbed AGN. This in turn will allow us, in particular, 
to study the absorption/reflection medium properties of the 
sample. These good quality spectra also permit a detailed 
investigation of the intrinsic continuum slope including the 
cut-off if present. The broad-band coverage available is a 
powerful and unique tool to address all these issues, making 
possible a simultaneous measurement of the intrinsic con- 
tinuum (both photon index and high energy cut-off) and of 
the reflected continuum together with the Fe line. 

Note that a number of sources in the sample have data 
already published in the literature both in the X-ray band 
(Sazonov et al. 2005, Shinozaki et al. 2006) and in soft 
gamma-ray range (Beckmann et al. 2006, Molina et al. 
2006). However, most of these works did not provide a de- 
tailed analysis of the broad band spectra of the sources in 
the sample nor attempted to put constraints on the high 
energy cut-off and reflection bump. The only exception is 
ESO 103-G35 analysed by our team using Beppo-SAX and 
INTEGRAL data; here we use a recent XMM-Newton ob- 
servation of the source in conjunction with an updated 
INTEGRAL spectrum in order to have in the sample a 
well studied object to use as a source of reference. 

In Sect. [5] we present the data of our small sample while 
the models used to fit the spectra are described in Sect. [3] 
Results and discussion are presented in Sect.2]and draw our 
conclusions. The single cases are discussed in more details 
in the appendix. 

2. Observations and data reduction 

2.1. INTEGRAL data 

The AGN discussed here have been observed with IBIS 
(Ubertini et al. 2003), the imager on board INTEGRAL 
(Winkler et al. 2003). This coded mask instrument is made 
by the combination of two detector layers: ISGRI (Lebrun 
et al. 2003), an upper CdTe layer sensitive in the range 

1 6 other type 2 AGN from Bassani et al. survey are be- 
low the assumed threshold flux: 3 (IGR J12415+5750, IGR 
J20286+2544 and IGR J17513-2011) had no X-ray available 
data when this analysis started while 3 (NGC 1068, Mkn 3 and 
NGC 6300) had BeppoSAX observational results published but 
all are peculiar objects the first two being Compton thick AGN 
and the third one the prototype of "changing look" Seyfert 2 ( 
i.e. a source that goes from thick to thin and viceversa). 



between 15 keV and 1 MeV, and PICsIT (Di Cocco et al. 
2003) a bottom Csl layer sensitive in the range 200 keV to 
8 MeV. In the present paper, we refer to data collected by 
the first layer only, since the sources are too weak above ~ 
100 keV for detection by PICsIT. 

For each source, INTEGRAL data from several point- 
ings performed between revolution 12 and 429 were added 
together in the following way to provide sufficient statis- 
tics. First ISGRI images for each available pointing were 
generated in 13 energy bands using the ISDC offline scien- 
tific analysis software OSA version 5.1. Then, count rates at 
each source position were extracted from individual images 
in order to provide light curves in the various energy bands 
sampled; since the light curves did not show any sign of 
variability or flaring activity, average fluxes were then ex- 
tracted in each band and combined to produce the source 
spectrum. A detailed description of the source extraction 
criteria can be found in Bird et al. (2006, 2007). We have 
tested that this method of spectral extraction is reliable by 
comparing the Crab spectrum obtained in this way with 
the one extracted using the standard spectral analysis. 

In Table Q] we report the details of the IBIS/ISGRI 
observations; best fit positions, exposures and count rates 
in the range 20-100 keV. The positional uncertainty for 
sources of this intensity is around ^3 arcmin (Bird et al. 
2006). 



2.2. Chandra data 

Chandra observed LEDA170194 on July 25, 2005, and 
IGR J12026-5349 and IGR J07565-4139 on June 16, 2005, 
with ACIS-I. Details on the observations are given in Table 
[T] Data were reduced using CIAO v3.2, following stan- 
dard procedures. It is important to note here that although 
Sazonov et al. (2005) used Chandra spectra in conjunction 
with INTEGRAL data, they did not attempt to estimate 
cut-off energy and reflection bump in their analysis proba- 
bly due to the limited amount of spectral information avail- 
able above 10 keV (a single bin in the 17-60 keV band). All 
Chandra sources suffer from pile- up. Using WebPIMMSQ 
we estimated the fraction of pile-up to be around 36 per 
cent for IGR J07565-4139, 30 per cent for IGR J12026-5349 
and >20 per cent for LEDA170194. The effect of pile-up is 
clearly visible through an hardening of the spectrum above 
4.5-5 keV. Regarding the pile-up problem, we point out that 
we have used a different approach than the one adopted by 
Sazonov et al. (2005): these authors removed the central 
pixel, i.e. excluding 80% of the instrument PSF and conse- 
quently used data up to 5 keV; we have instead maintained 
the central pixel and dealed with the spectra affected by 
pile-up through the PILE-UP model available in XSPEC 
which can correct the total spectrum for the pile- up effects, 
provided that this is less than 70% (Davis 2001). 

2.3. XMM-Newton data 

The XMM-Newton EPIC observation of ESO 103-G35, 
IC 4518A and IGR J10404-4625 have been also. Details on 
the observations are given in Table [TJ Data reduction was 
performed with version 7.0 of the SAS software, employing 
the most updated calibration files available at the time of 

2 http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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the data reduction (2006 March). Patterns to 12 (4) were 
employed in the extraction of the MOS (pn) scientific prod- 
ucts. In the case of IC 4518A pn and MOS spectra for the 
source were extracted from a circular region of 30 arcsec ra- 
dius while background spectra were taken from source-free 
circular regions of 10 arcsec radius. In the case of ESO 103- 
G35 MOS (pn) spectra were extracted from a circular re- 
gion of 35 (60) arcsec radius around the source while back- 
ground spectra were taken from source- free circular regions 
of 10 (40) arcsec radius. Finally in the case of ICR J10404- 
4625 MOS (pn) spectra were extracted from a circular re- 
gion of 40 (90) arcsec radius around the source while back- 
ground spectra were extracted from source-free circular re- 
gions of 20 (40) arcsec radius. Spectra are rebinned to have 
at least 20 counts in each spectral channel. Spectral fit 
were performed simultaneously with MOS1, MOS2 and pn 
in 0.5-10 keV. Cross-calibration constant pn/MOSl and 
MOS2/MOS1 were left free to vary and always found in the 
range 0.97-1.04. XMM-Newton data of ESO 103-G35 have 
been already published in Shinozaki et al. (2006). Their 
analysis was restricted to 0.5-10 keV energy range and is 
consequently not conclusive about Compton reflection and 
high energy cut-off. In Appendix [X] we will discuss also 
our results about absorption in comparison with those pre- 
sented in Shinozaki et al. (2006). 

2.4. ASCA data 

We finally analysed an archival ASCA spectrum of 
NGC 788. The relevant spectra and associated files were 
downloaded from the TARTARUS archivtS We fitted the 
spectra of all instruments SIS0, SIS1, GIS2 and GIS3 simul- 
taneously in the 0.6-10 keV band. The cross-calibration 
constant between the different instruments were left free 
during the fit and found in the range SISl/SIS0=1.06t^ 1; 
GIS2/SIS0=0.96t°;J°, GIS3/SIS0=1.09±g;Jjj. 

3. Spectral fitting 

The spectral analysis has been performed with XSPEC 
vll.3. Errors correspond to the 90 per cent confidence 
level for one interesting parameter (A% 2 =2.7). We fit- 
ted simultaneously the soft and hard X-ray spectra avail- 
able with INTEGRAL and either Chandra, or XMM- 
Newton or ASCA. The cross-calibration constant be- 
tween INTEGRAL and the low energies instruments 
{Cintegral) was fixed to 1. Possible miscalibration be- 
tween INTEGRAL and XMM-Newton, Chandra or ASCA 
(i.e. Cintegral 7^ 1) can mimic or hide the presence of 
a Compton reflection component above 10 keV. We choose 
to put this constant equal to one for all spectra because the 
fit we found with this value is good (all reduced \ 2 around 
1). However, in Sect. 2] we discuss the implications on our 
fit if Cintegral is left free to vary. 

To reproduce the broad-band spectra of our sample, 
we employed three different models. The first one (Model 
A) consists of an intrinsic absorbed power-law, a ther- 
mal soft X-ray component and a gaussian emission (Fe) 
line. The data and data/model A ratio for all source 
is shown in Figure [T] In the second model (Model B), 
we added to the primary emission a high energy cut- 
off; in this case the primary emission was modelled with 

3 http://tartarus.gsfc.nasa.gov 



an e-folded power-law (CUTOFFPL model in XSPEC). 
In the third model (Model C), we added to the previ- 
ous continuum a Compton reflection component; in this 
model the primary emission is modelled with an e-folded 
power-law plus reflection from an infinite cold slab. In 
the Compton reflection component (PEXRAV model in 
XSPEC, Magdziarz & Zdziars ki 1995D , the cosine of incli- 
nation angle of the reflector was fixed to 0.95. We assume el- 
emental abundances from Anders & Grevesse (1989). In the 
case of LEDA170194, IGR J07565-4139 and IGR J12026- 
5349 we choose to fix in all models the photon index r=1.7, 
i. e. the average value observed in AGN, the energy of the 
Fe line E^ e =6.4 keV and the intrinsic width of the line 
ape^lQ eV. These assumptions are needed because the low 
statistics of Chandra data did not allow us to constrain si- 
multaneously all the parameters of the models, T, E cut and 
reflection. 

The results of the fitting procedure with model A, B and 
C are presented in TableO Tableland Table[4]respectively. 



4. Results and Discussion 

All our models provide a good fit to all spectra, with % 2 /dof 
around 1. Model A, B and C have in fact an average reduced 
X 2 value of 0.97, 0.92 and 0.90 respectively. Here we present 
the main results of our analysis. 

4.1. The soft X-ray spectra 

All sources with the exception of IGR J07565-4139 show ev- 
idence for an unabsorbed component at energies less than 
2 keV. Recent high resolution spectroscopic data available 
with XMM- Newton on bright abso rbed AGN (Massaro et 
al. 2006, Guainazzi & Bianchi 2007) show evidence that the 
soft excess component is dominated by emission lines and 
originates in a gas photoionized by the primary continuum 
in which the density is decreasing with radius like r" 2 , as 
observed in the Narrow Line Region (NLR). On the other 
hand high resolution images available with Chandra show 
that in some cases (e.g. NGC 4549 Schurch et al. 2002) 
the soft X-ray emission has a thermal origin and is domi- 
nated by a starburst component. We have therefore taken 
two approaches, in the first we modelled the soft excess 
with a thermal component, i.e. a black-body model. With 
only one exception, the temperatures we found are in the 
range 0.2-0.9 keV (see Tabled]); in the case of IGR J12026- 
5349 the value of kT=1.7tQ2 keV is much higher than typ- 
ically observed in type 2 AGN. In the second, we investi- 
gated the alternative scenario in which the soft X-ray excess 
is reproduced by a scattered component, i.e. a power-law 
having the same photon index of the nuclear continuum 
but a different normalization. In this case we still obtain a 
good fit in the case of LEDA170194, NGC 788, ESO 103- 
G35, IC 4518A and IGR J10404-4625 (x 2 /dof are 25/25, 
291/298, 1081/1102, 102/126 and 1136/1071 respectively to 
be compared with those in Table [2]). The only exception is 
again IGR J12026-5349. The value of the ratio A IC /A soft 
(relative to the case where the soft and hard photon in- 
dices are frozen to the same value) for all sources is listed 
in Table [H This value provides an estimate of the rela- 
tive strength of the soft excess with respect to the primary 
emission. The ratios in Table [2] are of the order of few per 
cent for those sources that are well reproduced by a soft 
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Table 1. Journal of the observations 
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power-law component; this value is also in good agreement 
with that measured in a sample of absorbed AGN observed 
with XMM-Newton and Chandra (Bianchi et al. 2006). This 
soft component is likely produced by electron scattering of 
the intrinsic continuum and is observed in most Seyfert 2 
galaxies (Antonucci 1993, Matt et al. 1997). The evidence 
that the r 'soft and Thard are consistent with each other and 
that A/c/Aso/t is in the expected range of values, strongly 
support this scenario. 

In the case of IGR J12026-5349 a scattered component is 
not able to reproduce the soft X-ray excess (x 2 /dof=119/51 
when we substitute the thermal component with a power- 
law). Fitting with the scattered model Chandra data only, 
we get a good fit, however the ratio between the normal- 
ization of the two power-laws is ~ 0.12, making this sce- 
nario unlikely. The high value of kT and A/c/A so /i in 
IGR J 12026-5349 give indication that the contribution of 
the soft excess component to the emission is higher than in 
the other sources. This could be due to a starburst contri- 
bution. We will come back to this source in the appendix. 

If left free to vary, the photon index of the soft compo- 
nent remains consistent with that of the primary emission 
in the case of LEDA170194, ESO 103-G35 and NGC 788 
(r so/t = 1.018;?, r *oft = 1.8183 and r so/t = 2.1±8:i re- 
spectively to be compared with the values in Tabled]). This 
supports the "scattered" scenario for these three sources. 
On the contrary for IC 4518A if left free to vary the pho- 
ton index of the soft unabsorbed component is steeper with 
respect to that of the absorbed one (T so f t — 2.7l ; 4 ) sug- 
gesting that in this case the soft X-ray component could be 
of thermal origin. In the case of IGR J10404-4625 the soft 
photon index gets instead harder than that of the nuclear 
component (T so f t = 0.91q' 2 ); such a hard value is probably 
related to the presence of the Compton hump above 10 keV. 
This clearly shows the importance of fitting simultaneously 
the data from 0.5 to 200 keV. 

IGR J07565-4139 is a peculiar case, as in this source 
we do not find evidence for "continuum" soft X-ray emis- 
sion. A narrow emission line at 2.46±g;g| keV is required in 
the model for this source; the line has an equivalent width 
of 184l®g eV with respect to the unabsorbed continuum 
(see Table 2] model C represents the best fit in the case of 



IGR J07565-4139). The energy of the line suggests that it 
could be due to K a fluorescence of highly ionized elements 
like S or Si (SXV with the theoretical value of energy of 
2.430 keV and SiXIV with the theoretical value of energy of 
2.461 keV) and could be produced in a scattering warm gas 
photoionized by the primary emission as already observed 
in other obscured AGN such as NGC 1068 and Circinus 
galaxy (Guainazzi et al. 1999, Massaro et al. 2006). The 
lack of any other emission feature due to high ionized ele- 
ments prevents the creation of a ionization structure for this 
gas. If we assume that the warm medium has N# ~ 10 22 
cm~ 2 and ionization parameter Ux (defined as the ratio 
between the density of ionizing photons and the gas den- 
sity) in the range in 0.1-5, we can produce the observed 
line. This gas could be associated to the NLR and/or the 
[OIII] ionization cone s as demonstrated t o be the case for 
other obscured AGN (jBianchi et al. 2006|) . 

4.2. Absorption 

A column density in excess to the Galactic value was mea- 
sured in all sources, as expected in type 2 Seyfert galaxies. 
The minimum and maximum values of N# in our sam- 
ple are 0.1 and 40xl0 22 cm -2 respectively, with an average 
value N H =(10.6t°ol) x 10 22 cm" 2 (see Tabled}. For all our 
sources we have also estimated the Fx /F[OIII] A5007 ratios 
(see Table [5). The flux of [OIII] A5007 has been corrected 
for intrinsic galactic reddening following the prescription of 
Bassani et al. (1999) and is used here as an isotropic in- 
dicator of the intrinsic brightness of each source. We find 
that all Fx /[OIII] ratios are indeed well above those mea- 
sured in Compton thick sources and comparable with those 
observed in type 1 AGN (Bassani et al. 1999). The column 
densities measured as well as the Fx/F[OIII]A5007 ratios, 
shown in Table [51 allow us to confirm the already estab- 
lished Compton thin nature in the case of LEDA170194, 
IGR J07565-4139, IGR J12026-5349 and ESO 103-G35 and 
to establish that also NGC 788, IC 4518A and IGR J10404- 
4625 are Compton thin AGN. 

Malizia et al. (2007) recently analysed a sample of AGN 
observed with INTEGRAL and Swift/BAT providing a 
diagnostic tool (Flux (2-10 keV)/F(20-100 keV) softness 



Table 2. Model A composed by an absorbed power-law component, a thermal soft emission and a gaussian iron line. 





l N gal 




r 


A so ft/ AlC 


2 KT so /t 


2 E;ine 


2 __ 
CHine 


3 EW, illc 


r 2-10fceV 


4 F20-100fceV 


xVdof 


LEDA170194 


0.0369 


4.411:3 


1.7* 


0.01 


o-9±8:i 


6.4* 


0.01* 


< 1000 


2.7 


4.9 


16/24 


IGR J07565-4139 


0.458 


0.4tj>;» 


1.7* 






6.4* 
2.4±0.5 


0.01* 
<0.2 


< 548 
109l 7 g 


0.8 


1.5 


68/51 


IGR J12026-5349 


0.164 


qi+120 


1.7* 


0.12 


l-7l° -i 


6.4* 


0.01* 


< 3000 


10 


3.9 


57/50 


NGC 788 


0.0217 


401* 


1 6+ ' 1 


0.01 


0-301S-.S? 


6.3l° :l 


0.27t°;£ 


9 4 n + 177 


1.7 


3.8 


304/297 


ESO 103-G35 


0.0764 


18.81°:? 


1 8S+ ' 06 


0.002 


o-28±8:gS 


6.3l° :i 


0.4712:11 


19711 


6.1 


7.9 


1078/1101 


IC 4518A 


0.09 


15.01^ 


1 - J -0.1 


0.03 


n 99 +0.03 


6.261HI 


< 0.3 


qqc+197 
00<J -166 


0.58 


2.9 


76/125 


IGR J10404-4625 


0.137 


^•-■--o.i 


oq+0.06 


0.009 


n oq+0.08 


6.2ii° ; « 5 7 


19+ 011 


209±£ 


1.7 


2.5 


1132/1070 



(1) In 10 22 cm- 2 . (2) In keV. (3) In eV. (4) In 10" 11 erg cm- 2 s _1 . * Frozen during the fit. 



Table 3. Model B composed by an absorbed power-law component with high energy exponential cut-off, a thermal soft emission and a gaussian 
iron line. 
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Table 4. Model C composed by an absorbed e-folded power-law component, a reflection from a cold slab, a 
gaussian iron line. 
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Fig. 1. Spectra and data/model ratio for all sources of our sample. The ratios are plotted using our model A: absorbed 
power-law+thermal emission+iron line. 



ratio) to isolate peculiar objects and to investigate their able to confirm again the Compton thin nature of all our 
Compton thin/thick nature. Using this new tool they were sources, but identified IGR J07565-4139 as a peculiar ob- 
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Table 5. Fux of [OIII] A5007 corrected for absorption and 
its ratio with the unabsorbed flux in 2-10 keV. 
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IMasetti et al. 2006cl (3) Vace li et al. 1997. (4) Bassani et al. 
1999. (5) IMasetti et al. 2006bl 



ject as its location in this new diagnostic plot was outside 
the expected trend. 



4.3. The reflecting medium, the origin of the Compton bump 
and of the iron line 

We measure the reflection fraction through the pa- 
rameter R that represents the ratio between the re- 
flected and primary components. The model (PEXRAV, 
Magdziarz & Zdziarski 1995) assumes a perfect reflection 
from an infinite slab covering 2ir angle. This means that in 
the case of a face-on line of sight, R has to be in the range 
0-2. The importance of the reflection component is given 
by a comparison between Tableland TablelU The addition 
of one interesting parameter, R in model C, with respect 
to model B is statistically required in the case of ESO 103- 
G35, IC 4518A and IGR J10404-4625, with a significance 
of 99.98, 99.95 and 99.93 per cent respectively, according 
to the standard F-test. In the case of IGR J 12026-5349 and 
NGC 788, R is constrained but with a F-test probability 
of only 93.6 and 99.1 per cent respectively. In the case of 
LEDA170194 and IGR J07565-4139 we are only able to 
find an upper limit to R even if with an unphysical value 
greater than 2. We will return to this point later. In Figure 
[2] we plot the reflection fraction R versus the photon in- 
dex r. The values of R we obtain are in good agreement 
with those found in the Beppo-SAX sample of type 2 AGN 
(|Risaliti 2 002). Beppo-SAX average values are also shown 
in Figure H IC 4518A and NGC 788 show evidence of a 
flatter photon index than the other sources, this effect will 
be discussed in Sect. 14.41 and Appendix DO No evidence for 
the correlation claimed by Zdziarski et al (1999) between 
R and T is found (correlation coefficient r = 0.0582, sug- 
gesting no real evidences against the null hypothesis). It 
is important to stress here that the value of the Compton 
reflection fraction R is strongly dependent on the value of 
the cross-calibration constant Cintegral between the soft 
gamma (INTEGRAL) and X-ray data (Chandra, XMM- 
Newton and ASCA), that we assumed to be equal to one in 
all spectra. It is important to note that Kirsch and col- 
laborators (2005) analysed the Crab spectrum observed 
with different instruments, and they concluded that at 
20 keV the cross-calibration XMM-Newton/ INTEGRAL 
Chandra/ INTEGRAL and ASCA/ INTEGRAL was close to 
1 (within fews per cent). This strengthens the assumption 
we made on Cintegral- However we also checked our as- 
sumption a posteriori. In the case of XMM-Newton data, 
if left free to vary in the model C, the cross-calibration 
constant MOS1/ISGRI ranges between 0.5 and 2. In the 



Fig. 2. Reflection fraction versus photon index in our sam- 
ple. Dot-dash box show the avera ge values foun d in the 
Beppo-SAX sample of type 2 AGN (jRisaliti 2002j) . 



case of ASCA data of NGC 788, Cintegral is 0.8+^. In 
the case of Chandra data, Cintegral was much less con- 
strained with respect to XMM-Newton and ASCA, mainly 
because the larger gap between IBIS and ACIS (7-20 keV). 
In fact for Chandra data, if left free to vary in model C, 
the cross-calibration constant is found in the range 0.5-4; 
however the photon index fixed to 1.7 in the fit, still al- 
low us to measure the reflection fraction R in IGR J12026- 
5349, while only upper limits are found in LEDA170194 and 
IGR J07565-4139. In table H we report the value of R we 
find leaving Cintegral free to vary, and it is evident that, 
even if the value of R are less constrained, we are able to 
measure the Compton Reflection in four objects, and their 
best fit value are not very different from those obtained 
with Cintegral—^- We stress also that the errors associ- 
ated to R are higher than those of T or kT and this is due 
to the degeneracy between the different spectral parame- 
ters, in particular between T, E c and R. If we keep fixed T 
and/or E c , the values of R are much more constrained with 
respect to those reported in Table Q] and Fig. [2] 
The simplest scenario to explain the presence of the 
Compton hump above 10 keV in type 2 AGN, is one in 
which the reflecting and absorbing medium are the same, 
probably the putative molecular torus. Within this scenario 
the question to address is if the reflection fraction we find 
is consistent with the absorbing column densities we mea- 
sured. For N# = 10 23 — 10 24 — 10 25 cm~ 2 , the contribution 
of the torus to the 30 keV flux is 8, 29 and 55 per cent 
respectively (Ghiscllini et al. 1994). This means that the 
value of R measured in our sources is too high to be associ- 
ated only to the absorbing gas of the torus (see Table H|). A 
possible solution to this problem is that the absorber is not 
homogeneous and that a Compton thick medium covers a 
large fraction of the solid angle but not the line of sight 
as proposed by Risaliti (2002). This absorber can be iden- 
tified with the clumpy torus recently proposed to explain 
infrared observations of AGN and their classification in the 
unification theory scenario (Elitzur & Shlosman 2006). 

In the case of NGC 788, ESO 103-G35, IC 4518A and 
IGR J10404-4625, it has been possible to measure the Ka 
iron line. For the threee Chandra sources, LEDA170194, 
IGR J07565-4139 and IGR J12026-5349, we have freozen 
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Fig. 3. The iron line profile when the continuum is reproduced with model C. For NGC 788 stars, triangles, circles and 
squares represents SISO, SIS1, GIS2 and GIS3 respectively For IC 4518A, ESO 103-G35 and IGR J10404-4625 triangles, 
circles and squares represent MOS1, MOS2 and pn data respectively. 
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Fig. 4. Contour plots 99, 90, 68 per cent of the intrinsic width a vs Energy for the iron lines. 
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the energy of the (narrow) line Ei? e =6.4 keV and we were 
able to find only upper limits on the equivalent width. For 
this reason we present our conclusions about Fe line on 
the basis of the results we found in the other four sources. 
The profiles of the line when the continuum is fitted with 
model C, are shown in Figure [3l The energy of the line for 
all the objects is consistent with neutral or little ionized 
iron. The line is narrow in all objects as shown by the con- 
tour plot of a vs Energy (see Fig. [4| obtained fitting data 
with Model A. In the case of ESO 103-G35 the contours 
plot are only marginally consistent with a narrow line pro- 
file (cr=0.47t£;" keV with Model A, see Table©, however 
when the continuum is well reproduced with Model C, in 
which all the spectral components are modelled, we find an 
upper limit to the intrinsic width of the line (a < 0.06 keV, 
see Table 0] a Figure [3|) consistent with a narrow profile, 
as already observed by Beppo-SAX (Akylas et al. 2001). 
This underlines the importance of a good modeling of the 
continuum to investigate the line properties. 

In Figure [5] we plot the equivalent width versus the col- 
umn density for the four sources where the Fe line has been 
measured. If the Fe line is produced far away form the cen- 
tral source, a higher value of the column density will ab- 
sorb a higher fraction of the continuum at the energy of 
the line, but not the line photons; the total effect will be 
a larger value of the equivalent width. The errors on EW 
are too large to draw solid conclusion about a possible cor- 
relation, however fitting the EW data shown in Figure [5] 
with a constant K we found a value K=285±57 eV, with 
X 2 /dof=1.2. The same fit in the region below 3xl0 23 cm -2 
gives K=268±60 eV, with x 2 /dof=l (this best fit constant 
it plotted in Figure [5] with a dashed line). The only data 
point with N H > 3xl0 23 cm" 2 is that of NGC 788 which 
has a N# significantly higher than the other three sources. 
Even if we do not have enough statistics to draw any firm 
conclusions, we stress that the trend we find confirms the 
results of Risaliti (2002) on a sample of Seyfert 2 observed 
by Beppo-SAX. Beppo-SAX Njj average values below and 
above 3 xlO 23 cm -2 are also plotted in Figure [5) This be- 
haviour represents strong evidence that the line has to be 
produced (at least in part) far away from the accretion disc, 
being absorbed in a different way than the continuum. The 
evidence of a narrow profile (see Fig. 0]) strongly supports 
this scenario, and confirms that the best candidate for the 
line production site is the absorbing gas, probably a non 
homogeneous obscuring torus. 

The question is now if the absorbing gas alone, the torus, 
is able to produce the observed iron lines. The value of 
Nji in the four sources where we measure the iron line is 
in the range (4.2 — 36) x 10 22 cm -2 with an average value 
(16.9+g; 2 ) x 10 22 cm" 2 (see Table gj). A torus with column 
densities in this range could produce a fluorescence iron 
line with an average equivalent width in the range 10-200 
eV and with an average value of ^100 eV (Ghisellini et al. 
1994). These values are well below those we measure sug- 
gesting that the line has to be produced, at least partly, 
in a gas that is different from the absorbing one. Another 
possibility is, again, a non homogeneous absorber not cov- 
ering the line of sight uniformly. The upper limits of the 
equivalent width we found in the three Chandra spectra 
also support this conclusion. 

The measured values of the reflection components in our 
sample (both iron line equivalent width and reflection con- 
tinuum) suggest that the absorption is more effective than 
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Fig. 5. Absorbing column density versus Fe equivalent 
width. The dotted line represent the best fit constant in 
the region N# <3xl0 23 cm -2 . The box regions represent 
the average values of N# (la errors) found in the Beppo- 
SAX sample of type 2 AGN (| Risaliti 2002[) . 



the reflection, e.g. making the hypothesis that the absorb- 
ing/reflecting medium is not uniform, likely a clumpy torus 
(Elitzur & Shlosman 2006). Another scenario could be that 
the source is observed throughout a torus with a very shal- 
low angle, in this case the flux directly seen by the observer 
is due to the inverse Compton interaction of the hot plasma 
electrons off the seeds photons. This high energy radiation 
exiting the system with a smaller angle with respect to the 
torus plan, will interact with a sort of " grazing incidence" 
with the upper border of the torus and will see a very high 
column density (almost infinite). Conversely, the Compton 
scattered/reflected photons escaping in the direction of the 
line of sight, will have a Compton cone with a moderate 
optical path in our lign of sight. This do not apply to the 
lower energies, fully absorbed anyway. 

4.4. The intrinsic continuum and the high energy cut-off 

We measure the high energy cut-off in five out of seven 
sources in the sample (see Table 0} , and we find a lower 
limit to E c in the case of LEDA170194 and IGR J10404- 
4625. All these values (either measured or lower limits) are 
within the range found in other Seyfert galaxies (Perola et 
al 2000). 

The phenomenological model we employ to fit the in- 
trinsic continuum (i.e. power-law with a high energy cut- 
off), provides a good description of a two-phase model 
involving a hot corona emitting medium-hard X-rays by 
Comptonization and a cold optically thick layer (the disk) 
that provides the soft photons for the Comptonization 
(Haardt & Maraschi 1991; Haardt, Maraschi & Ghisellini 
1997). In Fig. [6]we plot the high energy cut-off vs the pho- 
ton index. It is important to note that the value of E c to 
be considered are those listed in Table 0] because above 
10 keV the presence of Compton reflection can modify the 
continuum slope and, in turn, the value of the high energy 
cut-off (we stressed in the previous sections that there is 
a strong degeneracy between the reflection fraction R and 
E c ). Fig. [5] shows that all the values we constrain are below 
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Fig. 6. High energy cut-off in keV versus photon index 



300 keV, strongly suggesting the presence of a high energy 
cut-off in the spectra of Seyfert 2 galaxies. 

In NGC 788 and IC 4518A we find a value of the photon 
index lower with respect the average observed in Seyfert 
galaxies. The flat photon index can be either due a bad 
modeling of the primary continuum, and this will be dis- 
cussed for the two sources in Appendix [XJ or it could be 
real, as expected in the X-ray background synthesis model 
(Gilli et al. 2007|) . The latter hypothesis needs more statis- 
tics to be confirmed. 



5. Conclusions 

We have presented spectral broad-band analysis of a small 
sample of seven absorbed Seyfert 2 selected at E> 10 keV 
by INTEGRAL /IBIS. The combined 0.2-200 keV spec- 
tra allowed us to perform a detailed study of the emis- 
sion/absorption properties and to test possible correlations 
between different parameters (e.g. R vs T or E c vs T); in 
particular, for the first time in these sources a measurement 
of the reflection components (both Compton bump and iron 
line) and high energy cut-off has been performed. An a pos- 
teriori check of our analysis procedure, which is affected by 
the limitation of using non simultaneous X-ray soft gamma- 
ray data, is performed providing for one well known source, 
ESO 103-G35, a direct comparison between our study and 
previous ones with Beppo-SAX data. The values of T, R and 
Nh we found for this source are completely consistent with 
those obtained by Akylas et al. (2001) analysing simulta- 
neous broad-band Beppo-SAX observations. This evidence 
strengthens the validity of our fitting procedure and makes 
reliable our results. The main conclusions or our analysis 
are the followings. 

— The average value of the absorbing column density is 
Nh = (10.6to'2) x 10 22 cm -2 , and the range observed 
suggests a Compton thin nature for all objects in the 
sample, that is a new results in the case of NGC 788, 
IC 4518A and IGR J10404-4625. This evidence is also 
confirmed by the ratio F x /F[OlU\ A5007. 

— The continuum was reproduced with an e-folded power- 
law having a photon index in the range 1.3-1.9 and a 
high energy cut-offs well below 300 keV; this finding 
suggests that a cut-off is a common factor in Seyfert 2 



galaxies with values in good agreement with those found 
in type 1 AGN. 

— We measure the Compton reflection component in five 
sources (IGR J12026-5349, NGC 788, ESO 103-G35, 
IC 4518A and IGR J10404-4625) and found an up- 
per limit in two others (LEDA170194 and IGR J07565- 
4139). No evidence for a correlation between R and F 
is found. We also observed the iron line (in the case 
of Chandra spectra we found upper limits only on the 
equivalent width). Both reflection components are not 
immediately consistent with production in the cold ab- 
sorbing gas identified in the molecular torus. Both the 
reflection fraction R and the equivalent width of the line 
are too high to be produced in the gas with the observed 
column density. A possible solution is that the absorp- 
tion is more effective than the reflection, e.g. making 
the hypothesis that the absorbing/reflecting medium is 
not uniform, like a clumpy torus, or that the source is 
observed through a torus with a very shallow opening 
angle. 

— A soft excess component emerging below 2-3 keV 
is found in all sources with the only exception of 
IGR J07565-4139. This component is well reproduced 
by a thermal black body model having temperature in 
the range 0.2-0.9 keV, as typically observed in Sy 1. 
Alternatively a power-law with photon index equal to 
that of the illuminating continuum is also able to repro- 
duce this excess is some cases; the ratio of soft versus 
primary continuum is of the order of a few per cent, 
as typically observed in Sy 2. A different case is that 
of IGR J12026-5349 which has, in the thermal model, 
a temperature kT=1.7^o 2 keV and, in the scattered 
scenario, Ajc/A so / t =0.12; both value are higher than 
observed in the other sources and, typically, in Seyfert. 
This strongly suggests that the source is characterized 
by a larger soft excess (possibly due to a starburst con- 
tribution) that is not present in the other sources. 

— IGR J07565-4139 is the only source in the sample that 
does not show any evidence for continuum excess at 
low energies. A narrow emission line around 2.5 keV is 
instead found in the Chandra data. 
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Appendix A: Analysis of single sources 

We describe here the details of the analysis of each individual 
source, we also compare our results to previously published 
measurements. 



scattered component) could suggest the presence of a starburst 
component. Our broad-band analysis well constrains the high 
energy cut-off E c = 170^gQ° keV, and we measured also a relative 
reflection R=1.31q 4. 

NGC 788 The photon index measured with the model C 
is flatter with respect to the average value observed in Seyfert 
galaxies, this could be either intrinsic or due to the presence of a 
high energy cut-off below 100 keV. In fact when fitting the data 
with model A, we found a photon index F=1.6, fully consistent 
with that expected in type 2 AGN. Another possibility could 
be that the absorber is more complex than that proposed here, 
e.g. a partial covering. We were able to constraint the Compton 
reflection fraction R=0.9lJ'7- 

ESO 103-G35 INTEGRAL data, oi ESO 103-G35 have been 
analysed by Molina et. al (2007) jointly with Beppo-SAX data 
from an observation performed in October 1996. Their data, 
however, referred to the 2th INTEGRAL catalogue (Bird et al. 
2006), and had only a net exposure of 41 ks. Our results, which 
are obtained with a longer exposure, are in very good agree- 
ment with previous results. In particular the value they found 
for the high energy cut-off (68I25 keV), perfectly matches our 
value (see Table [3] and [3} . They found an upper limit to the 
reflection fraction R<1.9, again consistent with our result (see 
Table [D >. XMM-Newton data of ESO 103-G35 have been anal- 
ysed bv lShinozaki et al. 20061 in 0.5-10 keV energy range. With 
their limited energy band it was not possible to obtain informa- 
tion about cut-off and reflection component, however the value 
of the intrinsic column density and the continuum shape they 
found is in very good agreement with the results of our broad- 
band analysis. We stress that the possibility to cover the energy 
range above 10 keV allowed us to model all the spectral compo- 
nents and then to extract the intrinsic continuum shape. 

IC 45 18 A In this object the presence of the reflection 
component above 10 keV is statistically required by the data. 
However the best fit value larger than 1 suggests that the ge- 
ometry of the reflector has to be more complex than that used 
in this work. The large contribution of the reflection component 
above 10 keV is also confirmed by the fact that the fit with 
model B, without reflection component, gives a very flat photon 
index. 

IGR J 10404-4625 The presence of the reflection compo- 
nent above 10 keV is statistically required by the data, with 
R=l.lio'.7- We were able to find a lower limit only to the high 
energy cut-off, E c >120 keV. 



LEDA170194 The fraction of pile-up estimated with 
PIMMs was around >20 per cent. This absorption column den- 
sity we found with our broad-band analysis is in good agreement 
with the result presented in Sazonov et al. (2005). In this source 
the presence of a high energy cut-off and a reflection component 
is not required (see Table [2] [3] and |4]) . We found a lower limit to 
E c of ~210 keV and an upper limit to R of ~2. 

IGR J07565-4139 The fraction of pile-up estimated with 
PIMMs was around 36 per cent. The absorption column density 
found in this source with our broad-band analysis is fully consis- 
tent with the value found by Sazonov et al. (2005) in their anal- 
ysis. Our broad-band analysis shows evidence that this source 
is peculiar in many ways, as already described in Sect. 4.1. We 
measured the value of the high energy cut-off which is 40±^° 
keV, the lowest value of our small sample. This very low value 
is probably due to the high upper limit to the reflection compo- 
nent <7; this high value suggests some different geometry than 
that used here to model the continuum reprocessed component. 

IGR J12026-5349 The fraction of pile-up estimated with 
PIMMs was around 30 per cent. Our analysis provides the ab- 
sorption column density completely consistent with that pub- 
lished by Sazonov et al. (2005). The enhancement observed in 
the soft X-ray (either modelled with a thermal component or a 



